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Abstract 

The full quantum mechanical collapse of a small relativistic dust shell is studied analytically, 
asymptotically and numerically starting from the exact finite dimensional classical reduced Hamil- 
tonian recently derived by Hajfcek and Kuchaf . The formulation of the quantum mechanics en- 
counters two problems. The first is the multivalued nature of the Hamiltonian and the second is 
the construction of an appropriate self adjoint momentum operator in the space of the shell motion 
which is confined to a half line. The first problem is solved by identifying and neglecting orbits 
of small action in order to obtain a single valued Hamiltonian. The second problem is solved by 
introducing an appropriate lapse function. The resulting quantum mechanics is then studied by 
means of analytical and numerical techniques. We find that the region of total collapse has very 
small probability. We also find that the solution concentrates around the classical Schwarzschild 
radius. The present work obtains from first principles a quantum mechanics for the shell and 
provides numerical solutions, whose behavior is explained by a detailed WKB analysis for a wide 
class of collapsing shells. 
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I. INTRODUCTION 



It is sometimes believed that the study of symmetry-reduced models, or mini-superspace 
models are of some importance in the understanding of the Physics in the strong curvature 
limit. This would be true, for instance, in the Big Bang or at the singularity formed by 
collapse of matter. There are a number of reasons why one would like to study the problem 
of the quantum collapse of a localized accumulation of matter. For large accumulations 
one might wish to understand the complete evolution of a collapsing system. A star, for 
example, might be expected to collapse through its horizon to form a black hole. The usual 
argument is that this ends all possibility of knowing what happened to the original star. 
However, the process of Hawking radiation will eventually bleed energy from the black hole, 
and, if quantum gravity is not taken into account, the original hole will disappear, with the 
consequent loss of coherence of the original system. 

One suggestion that avoids this problem is to assume that any curvature singularities 
inside the horizon are artefacts of assuming that gravity is a purely classical system, and that 
quantum gravity could replace these by some quantum state that shows no such singularity. 
When Hawking radiation has reduced the black hole to the point where quantum gravity 
comes into play, one might expect to find some sort of quantum remnant that might solve 
the problem of evaporation. This problem, even if we simplify the system drastically before 
quantization, is of daunting complexity, involving successively, classical collapse, quantum 
field theory on a curved background, and finally full quantum gravity coupled to quantum 
matter. 

A simpler problem is to consider the pure quantum collapse of phenomenological matter 
uncoupled to any field that allows radiation. While in general this involves turning off 
natural couplings, there is one reduced system that allows us to avoid radiative modes, at 
least classically. This is the classical collapse of a spherically symmetric accumulation of 
dust. Birkhoff's theorem shows that there is no gravitational radiation, and the dust matter 
has no radiative component (even though one might expect the expulsion of shells of matter). 
An even more drastically reduced case is the collapse of a spherical shell of dust, where the 
matter is constrained to lie on a (dynamical) spherical shell, thus avoiding blow off problems. 
This system provides perhaps the simplest quantum collapse scenario, where the effects of 
quantum gravity are uncontaminated by possible radiative effects. In principle one could 
use this model to study the collapse of a large shell that would pass through its horizon 
while still classical and evolve into some sort of quantum object near the classical curvature 
singularity. Without Hawking radiation this quantum object would always be hidden from 
outside observers, but it might give us some insight into the avoidance of singularities due 
to quantum effects. 

Note, however, that if one treats the Schwarzschild mass of the shell as a c-number (as we 
will below) 0], the mass will probably have a limit beyond which there is no valid solution 
to the quantum problem. This phenomenon would be the analogue of the limit of the charge 
of the nucleus in the relativistic hydrogen atom problem It is possible to make rough 
estimates of this limit. Thus, if we consider a particle of mass m orbiting in the gravitational 
field of a shell of mass M, the particle will be quantum mechanical if the radius of the shell 
is less that the gravitational Bohr radius of the particle, = / GMm? . It also seems 
reasonable to demand that be larger than the Schwarzschild radius of the shell for the 
particle to be observable (note that this excludes the large-mass collapse scenario defined 
above), i.e. > GMjc? (we ignore factors of order one). It is usual to argue that single- 
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particle quantum mechanics is no longer valid if Oq is less that the Compton wavelength 
of the particle, h/mc. The second and third of these conditions both give V Mm / Mp < 1, 
where Mp is the Planck mass. For stellar mass shells and elementary particles this bound is 
strongly violated, and elementary particles are produced well inside the horizon. For shell 
of roughly asteroid size (10^'' gm), the bound is saturated when the shell radius is roughly 
the Compton wavelength of the electron. We would expect production of light particles well 
outside the Planck length, and our single-particle quantum mechanics would break down. 
However, since we exclude all fields, we might hope that our formalism would hold until the 
mass of the particles produced was equal to that of the shell. In that case, the bound would 
be saturated when the shell mass was equal to the Planck mass (and the Schwarzschild radius 
equal to the Planck length). Thus we expect, in order for quantum mechanics to make sense, 
a limit on the mass of the shell of roughly the Planck mass, similar to a sufficient condition 
for self adjointness of the quantum operators given by Hajicek |0 for one of the models we 
will discuss below. 

Keeping in mind these caveats, we will investigate the quantum mechanical collapse of 
a roughly Planck mass shell of dust, from a rather novel point of view of assymptotics. 
This problem has been considered by several authors @, ^, |^, ^ |^. In all of these articles 
the Israel junction conditions were used to obtain a classical equation of motion for the 
curvature radius of the collapsing shell relative to the proper time of an observer on the shell. 
Once the equations were obtained, several possible Hamiltonians leading to these equations 
of motion were proposed and used to study the corresponding quantum mechanics. These 
Hamiltonians are both local and non-local. A fairly complete study of a local Hamiltonian 
is reported in The main conclusion of the above cited works is the existence of a discrete 
spectrum of bound states and a set of continuum eigenstates. It was also shown that the 
ground state hovers away from the classical Schwarzschild radius [Q, as might be expected 
from our qualitative arguments. For the non-local formulation similar results have been 
obtained, but they are less complete 0, §|. In this case the spectrum again turns out to be 
mixed but there is no detailed analysis given of the eigenfunctions. All these quantizations 
have qualitatively similar spectral properties. 

It is also important to remark, however, that all the above work, excepting , is not based 
on the Hamiltonian arising from the spherical reduction of the action of general relativity. 
Thus the question of how to perform and study the canonical quantization remains open. 
An important step in this direction is provided by the work of Kuchar and Haji'cek P, 0. 
In this paper the authors obtain the classical equations of motion of the shell as the exact 
finite dimensional reduction of the full relativistic Hamiltonian. 

The purpose of the present work is to continue the program outlined in and investigate 
the quantum mechanics of the Hamiltonian of Haji'cek, Bicak and Kuchar. (A variational 
procedure for deriving such a Hamiltonian has also been considered recently in |^). The 
first difficulty one must overcome arises from the multivalued nature of the canonical Hamil- 
tonian. This leads us to adopt a procedure of quantization which takes into account the 
relative size of the action on the different sheets of the Hamiltonian. When only the dom- 
inant action is taken into account, the problem can be formulated in terms of a non-local 
Hamiltonian. In this formulation the factor ordering of operators imposed by the require- 
ment of dilational symmetry of the half line (which is the space of values of the shell radius) 
forces a definite choice of time, which is related to but not the same as the proper time on 
the shell. 

Our procedure results in a non-local Schrodinger equation with vanishing coefficients 
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at the classical "Schwarzscliild radius" of the shell. Since our quantum mechanics relates 
to the internal parameters of the shell, with no reference to the external geometry, we 
will use the term "Schwarzschild radius" to mean the point where the curvature radius of 
the shell is equal to a constant of motion that is the Schwarzschild mass of the classical 
exterior geometry. The wave function that solves our non-local Schrodinger equation can 
be interpreted as the probability density of finding the shell with a given curvature radius. 



It is interesting to observe that non-local operators often appear in wave problems []TU 



while the vanishing of the main operator at the "Schwarzschild radius" is analogous to 
critical layer absorption in fluid mechanics [Tl], so the equations considered here combine 



two features which in classical fluid problems are separate. For these reasons, modern ideas 
of asymptotics |]10|, |l^ can be used successfully to obtain accurate information about the 
dynamics generated by our new nonlocal Schrodinger equation. 

The analysis of the solutions of the above mentioned quantum equation for the shell radius 



is performed using a generalization to non-local problems of WKB ideas |[TT[. The use of 
WKB ideas in a similar context in General Relativity has been explored in |13|, |T5|, |16|, |17 . 
In [jl5[ expected values of interesting physical parameters are calculated using the WKB 
ground state. While in [jl3|, |16|, tunelling amplitudes are computed for simple non-local 
operators. In the present work we perform a uniform WKB analysis of the continuum and 
discrete spectrum in order to study the full initial value problem for a collapsing shell. 
We find that there is indeed a phenomenon analogous to critical layer absorption, which 
decreases the amplitude of the eigenfunctions and the scattering states beyond the classical 
Schwarzschild radius. Moreover, we give a complete qualitative description of the dynamics. 
The analysis is complemented by a detailed numerical study (based on classical ideas in wave 
propagation |T2[ adapted to the present situation) for a representative range of parameters. 
The numerics give solutions that closely match the dynamics that was predicted qualitatively 
by the asymptotic analysis, and the results show a decomposition of the initial condition into 
bound states (which concentrate outside the "Schwarzschild radius" ) and a small probability 
current leaking to infinity. We also show how the present model explains, due to the presence 
of a critical-layer like behavior, the concentration of the solution of an initial value problem, 
with initial values outside the "Schwarzschild radius" , in the region outside this radius. 

To complete our analysis, we compare our numerical solutions with those obtained for 
a simpler non-local ad hoc model p|, § and contrast their differences. Finally we give 
arguments for how the present formulation could be used to obtain local information (around 
the space-time trajectory of the shell) about the space-time geometry. 

The paper is organized as follows: In the second section we review the variational for- 
mulation of the classical relativistic problem. The third section studies the classical motion 
and solves the problem of canonically quantizing the multivalued Hamiltonian. This section 
ends with a canonical quantum formulation in terms of a non-local equation for the shell 
problem. In the fourth section the Hamiltonian operator is analyzed both asymptotically 
and numerically. The fifth section is devoted to the detailed numerical study of the dy- 
namics for representative values of the parameters and initial conditions. The last section 
is dedicated to conclusions and suggestions for further work. The non-standard details of 
matching for the WKB solution are sketched in the Appendix. 
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II. FORMULATION OF THE PROBLEM 



In order to make our presentation as self contained as possible, as well as to fix notation 
and conventions, we recall briefly some of the details of the Kuchaf-Hajicek 0, |^ construction 
of the canonical dynamics of gravitational shells. The geometry of the shell is described by 

ds^ = -K^de + R\t)d^l\ (1) 

where the function R{t) is the (area) radius of the shell and ^ = A with r denoting the 
proper time on the shell. 

The Lagrangian for the dust matter is taken to be of the form 

= j ^M{A-'f^ - A)dt, (2) 

where M = Ai^R^pit) is the total mass of the shell and T = —Ut is the temporal component 
of the dust velocity . 

The gravitational fleld Lagrangian is expressed in terms of the variables ^ = x ^ 
proper velocity - and F± = 1 — where the Schwarzschild mass M± can be different on 
each side of the shell. For the special case of a collapsing shell M_ = 0, so the space inside 
the shell is flat. Thus 

£g = I (M+t+ + [L])dt. (3) 

Here T+ is the Schwarzschild time which determines the embedding of the shell in the space 
to the right of the shell. The function L is given by 

L = AR[{F + V^)^-VA{\F\~^V)], (4) 

where 

A{r]) = sinh~^ r], for F > (5) 

= sgn(r7)| cosh~^ Ir^ll for F < 0. (6) 

The jump [L] is obtained by using F+ = 1 — ^ on the right of the shell and F_ = 1 on the 
left. 

The total action is given by 

C = C^ + C^. (7) 

It is now possible to obtain the equations of motion. Variation with respect to M gives 
A = T. Variation with respect to T gives M = 0. Variation with respect to T+ gives M_|_ = 
constant. Finally, variation with respect to M+ gives T+ = KF^^^/F"^ + V^. The equation 
of motion is obtained by varying with respect to R. The energy condition which determines 
the trajectories for a given proper mass M. of the shell is 



f-f.--v^^-^. (S) 

The relative signs of the square roots are obtained from the requirement that the shell 
motion generate a cup space-time [0]. 
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This equation is taken as the starting point in the treatments in Refs. @, ^ • In these works 
there are several possible Hamiltonians that produce the equations of motion (^. However, 
The resulting quantizations, as we have pointed out in the Introduction, are not obtained 
from an exact dimensional reduction of the full general relativistic Hamiltonian of the theory. 
This leaves open the question (since the quantizations are not unitarily equivalent) of which 
canonical quantization one should use, and how the properties of these quantizations differ 
from those of ad hoc procedures. In the first place, one must obtain the classical Hamiltonian 
from the Lagrangian (^). This is done in ||^ by finding the appropriate canonical momentum 
and then constructing, by a finite dimensional reduction, the Hamiltonian. 

The Legendre transformation is given in H. Note that what was done there is equivalent 
to the full ADM procedure. The gravitational action leads to a Hamiltonian constraint 
instead of a Hamiltonian. By choosing an internal time, we can solve the constraint for 
the momentum conjugate to our time, which gives us a straightforward Hamiltonian (the 
ADM Hamiltonian). Classically, of course, this procedure is exactly equivalent to solving 
the unconstrained system and applying the constraint. Quantum mechanically there is 
an enormous difference, one procedure leading to an ADM Hamiltonian and the other to a 
Wheeler-DeWitt equation. In the ADM procedure, there is no guarantee that the evolutions 
for different times and different Hamiltonians are unitarily equivalent. The Wheeler-DeWitt 
route leads to the well-known "problem of time" in the solutions and difficulties in defining 
probabilities. Since we will use the ADM procedure and two different times in this article, 
we will briefly discuss problems associated with our approach. For the moment we will just 
give the results from necessary for our purposes. 

First, the expression for the canonical momentum P = |^ must be inverted to obtain V 

as a function of P. It turns out that for < i? < 2M the inverse function is double valued 
(see Equations (C33)and(C34) in 0). This poses a question of interpretation when using 
P instead of V. In Ref. the total Hamiltonian constraint for the shell plus gravity takes 
the form 

H^ + M = 0, (9) 

where the gravitational super- Hamiltonian H'^ is written in terms of P and R. Since is 
the conjugate variable to r, the proper (or curvature) time on the shell, the natural time 
choice is r, which gives a Hamiltonian Ai = H'^ . In Ref. there are several possible 
values for H'^ depending on whether R is greater than or less than 2M. We must construct 
the appropriate single valued function from these for a viable classical mechanics. To do 
this it is only necessary to construct the two-sheeted surface joined across the branch cut 
< i? < 2M at P = 0. 

The patching condition across the cut (analytic continuation) is found by taking into 
account the fact that V does not change sign across the cut Thus we have for R > 2M that 
the Hamiltonian is given by 

1/2 




i/f = -v/2P|l-^-^/l-^cosh^l . (10) 



For < P < 2M and P > the Hamiltonian is given by 

(: 



1/2 
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FIG. 1: R-P plane with cut in < i? < 2M. 
where the argument is positive. This is joined smoothly at P = with 

where the corresponding argument is positive and P < 0. Likewise, the other connection 
goes from P < to P > through the other sheet. The points on the different sheets are 
then projected onto the R — P plane to obtain the classical orbit. This is illustrated in Fig. 

It is to be noted that the double valuedness does not stem from the square root but from 
the P ^ V transformation. Also note that in the Legendre transformation the signs of V 
and P are not necessarily the same. The velocity V does not necessarily change sign on a 
given orbit when P does (although in some cases it does), due to the transformation. Thus 
the Legendre transformation is only possible in terms of the double valued Hamiltonian just 
described. As we shall see in the following Section, this will have non-standard implications 
for the quantum theory based on this Hamiltonian. 

Notice that it is not necessary to interpret the geometrical variable r as time. The 
constraint (P) is simply a linear combination of the geometrical variables H"^ and M. . 



III. CLASSICAL ORBITS, ACTION AND QUANTIZATION 

In order to quantize we must first study the classical mechanics of the problem and provide 
a suitable formulation for the multiple valued Hamiltonian. For the classical mechanics we 
have to solve for the level lines H = E. We anticipate these level lines to lie on two sheets. 
More specifically, we take R > 2M. Thus = E gives for the level lines the equation 

cosh^ = /(P,P), (13) 
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FIG. 2: {R,P) phase plane given by (|T|), (0), ([T|) and 

+ R^- MR) 



f{E,R) 



R^/WR - 2M 



(14) 



Provided f{E,R) > 1, Eq. (|T^) can be solved for P real. There are two solutions, one 
with P > and the other with P < 0. As i? — >^ oo there is a real solution only for E"^ < M^. 
For <i E^ < 2M^ the orbits do not extend to infinity; they have a turning point at 
R*{E) given by f{E,R*) = 1, with 



R*{E) 



ME^ + ^M^E^ + E^{E^ - M"^) E^{M + E) 



2(^2 - M2) 



2(^2 - M2) 



> 



(15) 



for E'^ > M^. On the other hand, for E^ > 2M'^ there are no orbits which represent classical 
trajectories. For these orbits, P oo as i? — ^ 2M. This part of the orbit is labeled Oi in 
Fig. |. 

Let us now consider the situation when R < 2M. When R approaches 2M from below, 
we have that the orbits are solutions of 



sinh^| = /(i?,i?). 



(16) 



The solutions match as \P\ oo since in that limit | sinh -^1 ~ cosh Thus the orbits can 
be continued, and if we pick the collapsing orbits, i.e. P < they are continued as solutions 
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of (|T6D up to P = 0. We find that P = at 



«.«(^).M!±11!, (17) 

obtained by solving f[E,R{E)] = 0. This orbit is labeled O2 in Fig. |^. Now, at this point 
the solution passes to the other sheet and continues as a solution of 

|sinh^| = -/(E,P), (18) 

with P > 0. This is labeled as O3 in Fig. 0. The mirror image of the above description 
applies to the situation for expanding orbits. Finally for the case < the situation is 
the same as before but now the orbits extend as incoming or outgoing orbits as P — >■ 00. 
These orbits are labeled in Fig. ^ as O4. This description reproduces the classical collapse 
as well as expanding motion with turning points. This is to be contrasted with the phase 
plane obtained in |]T^ where the local coordinate system used gives only the portion 



of the Hamiltonian. There the unbounded orbits represent collapse up to a point and the 
bounded ones represent expansion and ultimate collapse. There is a separatrix dividing the 
phase of these two types of motion. Also, in the quantum mechanics developed in the above 
cited paper there is an exponentially small probability of total collapse. In our analysis, in 
contrast, the orbits do not have a separatrix and they all show ultimate collapse. We will 
show later that the collapse turns out to be localized around R = 2M via a phenomenon 
of critical layer, the nature of which is unraveled by the Hamiltonian that we use. This 
phenomenon is not captured by the local expansion used in With this picture we 

turn now to the problem of quantization. 

We will choose our internal time to be the proper time on the shell, r, and E is its 
conjugate variable. We must now take H'~^ to be an operator on an appropriate Hilbert space, 
and P to be a multiplication operator while cosh P/P and its functions are to be realized as 
a pseudo differential operator on the Hilbert space [^. However, in order to carry out this 
canonical quantization we require a single-valued Hamiltonian, while, as already noted, the 
true Hamiltonian of the shell given in Equations ([TT|) and (12) is multi- valued. This situation 



is new on both counts, and because of this we need to introduce a proposal, which we believe 
to be new, to quantize the system. First we will consider possible definitions of the operator 
cosh P/P. Hajicek proposes quantizing the shell radius, which is restricted to the half 
line P > 0, by considering the Hilbert space of functions \1/ in L^[0, 00) with respect to the 
measure R~^dR. The observables are taken to be P and H = —iR-^. The factor ordering 

is taken to make the operator cosh P/P formally self-adjoint with P = H/P with respect to 
the measure R~^dR. Since the hyperbolic cosine contains only even powers of its argument, 
it is sufficient to take a factor ordering that makes (H/P)^ formally self adjoint 0]. In our 
case, due to the 1/P factor, it is easy to check that this ordering is not possible. This 
forces us, in order to preserve symmetries due to the dilational symmetry of the half line. 



to redefine the cosh P/P operator as simply cosh P. The operator | sinhP| = ycosh P — 1 
can be shown to be self adjoint by means of spectral methods. 

This quantization leads us to a new choice of time (lapse A) in order to re-scale the 
momentum. The freedom we have to choose A allows us to construct a rescaled Hamilto- 
nian and a canonical quantization scheme. Thus if we now take A = /i/P in Eq. (P, the 
Lagrangian (^ becomes 

L = ^{{F + V^f/^ -VF{\F\~^'^V)] (19) 
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Note that the resulting new canonical momentum obtained from 

does not contain the factor 1/R. Finally, Eq. (||) is exactly the same, except that now 
Y — RR ^\{iq\^ allows us to recover R from V . The Hamiltonian with this new choice of 

lapse is essentially the same as in Eqs. ([T0|)-(p!2D, except that P /R is replaced by P. With 
our new time, we can obtain (as in Ref. 0) coshP as a formally self-adjoint operator. 
The second difficulty we must address is the question of how to handle the multi- 



valuedness of the Hamiltonian. By analogy to the path integral formulation we consider 
the classical action 

R 

PdR 







along an orbit, say O3, O2, or Oi. Clearly the action along the dashed, O3, portion of 
the orbit is small compared to the action along the O2 and Oi portions in solid line. This 



occurs because R*{E) » R{E). An approach similar in spirit was used in [12] where the 
action was approximated in the complex plane and then the path integral was evaluated by a 
saddle point approximation. Here we approximate the real action and obtain an approximate 
operator equation. Also, since the classical Hamiltonian involves a square root, it is necessary 
to define first the square and then take the square root using spectral calculus. To this end 
we begin by defining H as the square of if*^ in the form 

H = 2{{R^ - MR) ~ R^/^\R-2M\^^'^ cosh P} R>2M, (21) 
H = 2{{R'^ - MR) - R^/'^\R-2M\^^^\smhP\} R<2M, (22) 

once the small contribution of the O3 portion of the orbit is neglected This form of H 
introduces classical turning points at R{E) which are not present in the original Hamiltonian. 
The implications of this fact will be discussed below. 



We now turn to the quantization of (|2T| - p2D . As we have already remarked, the operator 



coshP is transformed into a self-adjoint operator on functions \& relative to the measure 
dR/R. A unitary transformation \l/ turns coshP into a self-adjoint operator in 

L^(0, 00) relative to the measure dx. The functions must have all their even derivatives 
equal to zero at x = 0. We thus have 



H^{x) = 2|(x^ -Mx)^ 

-x^/^ \x - 2M|i/2_ J pQgj^ ^ e^'=^^(fc, t)dk[ , 

X > 2M, (23) 
H^{x) = 2{(x2 -Mx)^ 

-x3/2|x-2M|i/2— y |sinhfc| e'^''<il{k,t)dk\ , 
X < 2M, (24) 

where the Fourier transform ^' of is defined by 



if{k,T) = I e-^'=^^(x,r) dx. (25) 

QO 
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Now, since H will be shown to be positive, the spectral calculus gives the final equation in 
the form 

= Hy\, (26) 

and u{x,0) = \E'(a;) is the given initial condition. This procedure gives us a well-defined 
equation for the quantum evolution of the radius of the shell. If one prefers the Wheeler- 
DeWitt approach, one need only reinterpret ( p6D as the equation that result from applying 
the Hamiltonian constraint (|^) as a quantum equation, and think of r as simply a gravita- 
tional variable. The solutions will be the same. 

Before we solve the equation we need to make several remarks. In the first place, note that 
pseudo-differential equations like arise in classical wave propagation |TD[, where the non- 



locality of the main operator reflects the incompressibility of the fluid. The zero coefficient 
at X = 2M in front of the main operator also arises in the study of wave propagation 
in stratified fiuids in the presence of inviscid critical layers WM. However, in the case of 



critical layers the main operator is local. For these reasons we will use the same analytical 
and numerical methods developed for those problems, with appropriate modifications to 
take into account the presence of the vanishing coefficient multiplying the main nonlocal 
operator in order to solve the problem posed in (pB]). 

One final remark is that we will treat M as a c-number similar to the charge on the electron 
in the hydrogen atom problem. Actually there is no reason for this since M is merely a 
constant of motion, that is, a phase space function or g-number. A simple analogue might be 
a particle of mass m moving in two dimensions, where in one of the directions, z, the particle 
is free, and in the other, y, moves in a harmonic oscillator potential. The eigenfunctions 
have the form e^^^^ipniUiPz)-, where the ipn are harmonic oscillator eigenfunctions. We may 
consider an eigenstate of Pz and study the evolution of 4'{y,t) with constant, but there 

is no reason to do so. One can just as well study more general functions of the form 

■ pi 

I f{Pz)G^^''^^~^^*'4'{yjPz)G~^^"'^dpz- In our case, we have taken an eigenfunction of M, 6{M — 
Mo), but one should study more general situation. 

Finally, we compare the solution of Equation (^Gj) with another model studied in P, §] , 

M^]_r^ cosh k e'^^'uik, t)dk - ^. (27) 
dt 27ri-oo ^ ^ X ^ ^ 

Note that further other models considered in the literature have similar structures with 
different symbols for the nonlocal operator, or are similar to (p7|) but with functions of 
local operators. The cosh k terms appear in all models, since a proper time variable is used 
in their description. The coordinate system chosen to cross the classical horizon produces 
the nonlocal terms. The potential term is common to all the models. However, the exact 
reduction leading to (^) has the new feature of having a vanishing term multiplying the 
main operator. We will discuss the implications of the similarities and differences between 
the models in the following sections. 



IV. WKB SOLUTIONS, SPECTRA AND THE QUALITATIVE BEHAVIOR OF 
THE SOLUTIONS 



As X — 2M the term in front of the nonlocal operator in (p3D -(^J) produces high fre- 
quency spatial oscillations. Because of the non-locality we look for an approximate solution 
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in the WKB form of (|^). We take 

u = A(x)e''(")e-^^" (28) 

The phase 6 satisfies the eikonal equation (^) where P / R is replaced by 6' (since P/R is 
replaced by P in the new formulation), and the derivatives of A are assumed to be small 
relative to 6'. The eikonal equation and amplitude equation are readily obtained from the 



average Lagrangian |T0, ITl 



J J [i{u* drU - u drU*) + u* Hu] dxdt. (29) 
The eikonal equation takes the form 

|,i„h«'| = z£V2±^!_Z^. ,<2M (31) 

The first equation can be solved for real 6' (which represents an oscillatory solution as 
X oo) provided f{E, x) > 1. This occurs for x > provided E"^ < M^. On the other hand, 
when > real solutions are only possible if 2M < x < R*{E) = ^M^- Beyond R* 
the solution becomes complex, and this region, ^ E^ < 2M^, provides candidates for 
the point spectrum. On the other hand, the region < E^ < provides candidates for 
scattering states. Note that E < is not in the spectrum. In the region R{E) < x < 2M 
where R{E) = i(M + y/M^ + 2E^), Eq. has real solutions. Even though d' goes to 
infinity a.s x 2M from above or below, 6' matches, since cosh 6*' ~ | sinh^'| as \6'\ oo. 
For < x < R{E), 6' is again complex. This produces an exponentially decaying solution 
as X ^ 0. This is a consequence of the approximation we have used, which replaces a small 
action by zero in the region < x < R{E). Also note that when E^ < there are no 
solutions defined in the whole interval. This shows that H is indeed positive. 

Given the phases, the amplitudes are determined from the average Lagrangian by varying 
with respect to 9'. Using 

J [H{e')AA* - EAA*]dx, (32) 

we find that 

^ {\A\^Hi9')) = 0, (33) 

which gives = constant/if (6''), where in each region 6' is chosen according to Eqs. (|30D 
or (HD. 

We now give the explicit expressions for the scattering states and for the bound states. 
We begin by finding the phases. For R{E) < x < 2M we have 



^i(x) = / In w 
Jr{e) I \ 

, ie - Mi 



(e -Mi- Ey2Y ^ ^ 



^3/2|^_2M|V2 




(34) 
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while for 2M < x we have 



02 = 0i(2M) + r InJ . 

J2M I \ 

- 2M)i/2 




2M) 

(35) 



Note that when E'^ < M^, (x) is defined for all x > 2M. On the other hand, when 
E"^ > M^, 02{x) is real only up to a; = R*{E). The amplitude associated with 6i is given by 

= [(-^72 + x2 - Mx)2 + x^\x - 2M|]V4' R<x<2M (36) 
while that associated with 6'2 is 

" [(-^2/2 + x2 - Mx)2 - x3|x - 2M|]V4- 2M < X (37) 

Note that when x = 2M, Ai = A2. Also, when < E"^ < M^, A2 is always well defined for 
X > 2M and as a; ^ 00 it is given by A2{x) ~ On the other hand, when E"^ > M^, as 
X — >• R*{E), ^2 ^ 00. This indicates, as usual [1^, [l^, the presence of a transition region 
which has to be included in order to construct the exponentially decaying eigenfunctions 
of the point spectrum. Finally, as x — * R{E), the amplitude Ai is finite. Due to the 
nonanaliticity of | sinh6''| at 6*' = when x = R{E), the behavior for x < R{E) will be 
shown in the Appendix to be algebraically small. With these considerations we can now 
construct the WKB functions which represent both scattering states and bound states. 

We begin with the scattering states. From the Appendix we have that the solution for 
X ~ R{E), but less that R{E), matches to the solution 

Ai{x)smei{x). (38) 

In turn, this continues into 

A2{x)sme2{x), (39) 



where Ai, 6*1, A2, 62 are given by Eqs. (p4D-(p^). Thus the function Lp{x, E) defined by Eqs. 
(^) and (^) is a solution of the eigenvalue equation that is not square integrable and, as 
expected, represents a perfectly reflected probability wave impinging from infinity. 

The bound states occur for E'^ > M^. In this already noted above, A2 ^ 00 as 

X R*{E). This is typical turning point behavior, and the appropriate transition function 
will be shown in the Appendix to be the Airy function. Beyond R*{E) the phase is complex 
and the solution decays exponentially. The matching condition between the oscillating and 
the Airy function gives the usual Bohr- Sommerf eld condition for the eigenvalue E. It takes 
the form 

e2[R*iE)] = in + ^)n. (40) 

This provides a transcendental equation for E. It is clear from (|40D , since 62[R*{E)] 00 as 
E'^ — >■ from above, that there are infinitely many eigenvalues accumulating at from 
above. An illustration of the W.K,B. eigenfunctions is given in Fig. (3). Notice that the 
oscillations are concentrated at the "critical layer," x ~ 2M. It is also important to notice 
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that the size of the support of the oscillatory part of the eigenfunction which extends up to 
R*{E) increases very rapidly even for low eigenvalues (n = 3,4). 

The eigenvalues of H are the solutions to (|3y) and the continuous spectrum extends from 
zero to M. The analysis just presented shows also explicitly that the eigenf unctions are 
localized around R = 2M because of the critical layer behavior. Again, the existence of a 
point spectrum here is to be contrasted with the results of where the spectrum is only 



continuous. This effect was also found in other model Hamiltonians but no intrinsic 
explanation was possible. Finally, our spectrum and eigenf unctions must be contrasted with 
those of other models. In all models, including the present one, the spectrum is always 
mixed. However, previous models have an unbounded continuous spectrum, and their point 
spectrum has a ground state If. This situation just reflects the fact that previous 

models were based on a choice of an ad hoc Hamiltonian that reproduced the equation of 
motion. Note, nonetheless, that in principle the inverses of these ad hoc Hamiltonians are 
equally legitimate choices, as they also reproduce the classical equations of motion, but their 
spectrum will have a bounded continuous part that is below the point spectrum. Thus even 
if the quantum problems are not unitarily equivalent, they have the same spectral properties 
as the ones obtained from the finite dimensional reduced relativistic Hamiltonian. On the 
other hand, the form of the eigenfunctions in those models was found to be similar to the ones 
for the hydrogen atom ^J, while for the Hamiltonian considered here the eigenfunctions 
and the scattering states have a large number of oscillations concentrated around x = 2M. 

In the following section we consider the quantum dynamics generated by the Schrodinger 
equation ( p6D and compare it with the dynamics of the non-local Hamiltonian considered in 
0, ^, 1^ and which leads to 

1^ = ^ r cosh k e''^" u{k, t)dk - ^, (41) 

However, before closing this section we would like to make some remarks concerning the 
spectrum problem associated with (|4l|) which show that its eigenvalues and eigenfunctions 
behave as hydrogen-like. Thus taking u = vexp {—iEt), the eigenvalue equation becomes 

Ev = — coshfce^'^^ v(k,t)dk - ^. (42) 
2n J-oo X 

The WKB eigenfunctions take the form 

V = Ae*^(^) (43) 

and the phase 6 satisfies 

coshe' = E + -. (44) 

X 



Consequently for E > 1 there are always two real solutions for Eq. (|4^) , and the case E > 1 
is always part of the continuous spectrum with eigenfunctions perfectly reflected at x = 0. 
Alternatively, when E < 1 there are turning points for 

E + (^ = l, (45) 

X 

which gives x{E) = jyze)- quantization rule for the point spectrum is 

r^i^) ( d\ 1 

/ cosh"^ \E + -\dx = {n + -)ti. (46) 

«/ \ Ou J 
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This equation shows that since x{E) — > oo as ^ 1, there are infinitely many eigenvalues 
accumulating from below at E = 1. Also for E = the integral is smaller than ^. This 
implies that the spectrum is bounded from below. We thus see that the eigenvalues and 
eigenfunctions of ( ^2] ) behave indeed as hydrogen-like (see also P). With this information 
we shall next consider numerical solutions of (^) and (^) and provide an interpretation of 
the results. 



V. QUANTUM DYNAMICS AND NUMERICAL SOLUTIONS 

This section has three parts. In the first one, we focus on the general structure of 
the solutions of the quantum problem and show that, in certain regime of interest, the 
quantum dynamics can be approximated by a simpler Hamiltonian. In the second part, 
we describe the numerical scheme used in the calculation. In the last part, we describe 
the numerical solutions for different initial conditions and compare them with the WKB 
solutions previously studied. 



A. Spectral properties and analytical results for the dynamics 

Let us begin by first remarking on the basic structure of the solutions of the Schrodinger 
equations for the non-local Hamiltonians (^) and (|^). The solutions of the initial value 
problem take the form 

oo 

uix, t) = J2 Cn^n{x, K)e-^^"" + / c{E)^{x, E)e-*^" dE. (47) 

n=l 

The Cn and ciE) are just the projections of the initial condition on the discrete and con- 
tinuous eigenfunctions respectively. The second term in (^) decays in time as r~^/^. It 
is important to notice from equation ( |5^ ) that all the states in the continuous spectrum 
have high wavenumbers in the vicinity of the critical layer. It is therefore expected that the 
radiation produced will be of short wavelength in the vicinity of 2M and of long wavelength 
as X — i> oo. Since the eigenfunctions are strongly localized between R*{En) and R{En), the 
first term is localized. Since ip ^ as x ^ the dynamics will show small probability of 
total collapse. 

We have shown that the discrete spectrum is concentrated around E ~ M, while the 
continuous spectrum only contributes to the decaying part of the solution. Thus the solution 
is dominated (in particular for large time which is the relevant behavior for the physics of 
collapse) for values of the spectral variable close to M. We can take advantage of this 



fact in order to simplify considerably the numerical integration of (pG]). To achieve this 

simplification use H^^^ = \J [H^/'^Y. Using the spectral representation in terms of the 
eigenfunctions, we have 



H^/^^{x) = ^El Cni^n{En, x) + c(E)^(E, x)dE. (48) 

Now, since the spectrum is concentrated around E'^ ~ M^, we have, using the linear ap- 
proximation around M^, 

1 
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This gives the approximation for H^l'^ in the form 

Thus, due to the nature of the spectrum, the equation 

idr^ = H^/^^ (51) 

is well approximated by 

With a change of variables -ip = e"*^^^/^!; (x, r/2M) the equation to be solved numerically is 

idtv = Hv, (53) 

where t = r/2M. This reflects the combination of two facts : first that classically H^^^ and 
H have the same orbits with rescaled time, and second that since the discrete spectrum only 
runs over a small energy range, the solution to (|51|) is well approximated by the solution to 



B. Numerical scheme and numerical stability 

We now outline the numerical scheme used for the solutions of the equations. The equa- 
tion to be solved is (|53|). These equation was solved numerically using a pseudo-spectral 



method based on that of []T2[ for the Korteweg-de Vries equation. This method involves 
finite Fourier transforms in the x-direction, with the equation propagated forward in the 
t-direction by fourth-order Runge-Kutta integration. The finite Fourier transforms are eval- 
uated using fast Fourier transforms (FFT). 

The partial differential equation (^) is numerically solved in the interval —L/2 < x < 
L/2. Let us denote the Fourier coefficients of v by and let ^„ = 2nn/L. Then the Fourier 
integral operators in (^) can be evaluated as 

1 n=m 

- coshen^)ne-*«"^ (54) 

n=—m 
-I n=m 

- Yl |sinhen|^)ne-*«"^ (55) 

n=—m 

respectively, where 2m + 1 terms are taken in the finite Fourier series. 

Once the Fourier integrals have been evaluated, the partial differential equation (^) is 
then integrated forward in time in Fourier space. This integration is carried out in Fourier 
space due to its superior stability. Let us denote the right hand side of (pBD by Gv. Then if 
F denotes the finite Fourier transform, we have 

^^F{v) = -zF{Gv). (56) 



The ordinary differential equation (|56|) is integrated using fourth-order Runge-Kutta inte- 
gration to find F[v). The solution v is then determined by inversion using the FFT. The 
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pseudo-spectral method applied to the partial differential equation (|53|) has an inherent 
numerical stability problem which will now be discussed. 

The interface at x = 2M is a boundary layer in which the solution for v rapidly oscillates. 
This rapid oscillation results in the formation of high frequency modes in the finite Fourier 
series. Then, due to the equivalence between small x in physical space and large ^„ in Fourier 
space, these high frequency modes generate instabilities for x ~ which then propagate into 
the rest of the numerical domain. These high frequency modes are a major problem for ( |55[ ) 
as the hyperbolic functions in the Fourier integrals magnify their influence. To reduce these 
numerical instabilities, the numerical method as outlined above was changed in two ways. 
The first was that the governing equation (^) was changed to 



idtv = 2[{x^ - Mx)v 

r 1 1/4 1 ~1 

-^^/^(x - 2M)2 + — / cosh ke'''''v{k,t)dk\ 

J 2n J-oo J 



X 

X > 2M, (57) 



idtV = 2{(x2 - Mx)v 

r 1 1/4 1 r°° 1 

r^/^ {x - 2Mf + — / I sinh k\ e'^''v{k, t) dk \ 



—X 

X < 2M. (58) 

The modification of the coefficient of the Fourier integral terms in the equations was done to 
help resolve the boundary (critical) layer. The system ( p7D and ( |5BD was solved as previously 
explained in (pB]). In the numerical solutions, e was taken to be 0.1. The second modification 
was that a low pass filter was applied when evaluating the integrals in (|F7| ) and (^B|). This 
low pass filter is equivalent to setting cosh^„ and sinh^„ to for large It is interesting to 
note that the concentration of the spectrum at -E ~ M allowed us to control the instability 
by using the very accurate linear approximation given in (^3|). In other problems it may 
turn out that the spectrum is not concentrated in a small region, in which case one must find 
the numerical square root of the matrix FG. This process must be carried out in a manner 
consistent with the filtering, and this procedure is currently being studied. The choice of the 
number of modes used to obtain an accurate numerical solution was dictated by comparison 
of the long time numerical solution with the WKB solution as an initial condition with the 
same W.K.B solution. 



C. Numerical solutions 

We begin by considering the simpler problem of solving Eq. (^) with a Gaussian as an 
initial condition. In Fig. ^we present the result for the integration of (|1) with the Gaussian 
initial condition 

«(a;, 0) = e-(^-^'«)'/\ (59) 

for the parameter values A = 1.0 and Xq = 7.0, which was reflected as an odd function about 
X = to satisfy the boundary conditions at x = 0. This figure clearly shows the decomposi- 
tion of the initial condition into smooth normal modes which are localized away from x = 0. 
Notice that the modes are smooth and resemble the hydrogen atom eigenfunctions, as can 
be seen from equations (H3) and (BB) and the representation in terms of Laguerre functions 
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FIG. 3: Numerical solution of model equation (^) for initial condition (59) with xq 
A = 1.0. Here /? = 0.5. t = 0: ~- ; t = 20: ■ t = 50: ■ ■ t = 100: 
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To study the problem posed by Eqn. (|^) we begin by assessing the validity of the WKB 
analysis. To this end we use the WKB eigenfunction of Equations (|38D and ( pQ]) for n = 2. 
This eigenfunction is shown in Fig. ^. In Fig. ^ we show the evolution under the Schrodinger 
equations ( pTf ) and (^) of this WKB solution as an initial condition. It can be seen that a 
small amount of radiation leaks to infinity and that the WKB solution is distorted by high 
frequency waves. The high frequency oscillations are due to the non-smooth WKB initial 
condition and these high wavenumbers then spread under the influence of the non-local 
nature of the governing equations ( |F7| ) and Ignoring these high frequency waves, the 

full numerical solution compares well with the WKB approximation. This good comparison 
gives confidence in using WKB solutions to explore general initial value problems. 

Taking into account the above considerations, let us consider the initial value problem 
for ( ^7|) and (|58D with a Gaussian initial condition. We have two cases: xq > 2M and 



Xq < 2M. Let us consider first the case xq > 2M. In Fig. 5 we show the evolution of the 
initial condition at t = 5, at which time the solution has reached the steady state. For 
M = 5 this figure shows how the initial value is decomposed into the eigenfunction and a 
small amount of probability current that is propagated towards infinity. This result can be 
accounted for by the following WKB analysis. 

The discrete eigenfunctions have their oscillatory part concentrated around 2M. Thus the 
Gaussian projected onto these modes decomposes into the oscillations shown in Figs. ^ and 
1^. The projection onto the continuous spectrum produces the small current which spreads 
to infinity with decaying amplitude. The high frequency of the radiated waves stems from 
the short wave behavior displayed in the WKB solutions. 

The situation for the opposite case Xq < 2M is shown in Fig. ^ The chosen mass parameter 
was M = 10 and the initial value was taken to be Xq = 4. Here we have the opposite situation 
occurring, namely that the evolution is confined to the left of x = 2M, while for xq > 2M 
we found that the solution was concentrated to the right of x = 2M. This observation may 
be explained in terms of the WKB amplitude given in (^) and (|57D. From those equations 
it follows that the amplitude A2 is small at x = 2M and this, in turn, blocks the passage 
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FIG. 4: WKB solution given by (||) and (|9|) with M = 10 and n = 3. 
0.4 I , , , , 1 




FIG. 5: Evolution of the WKB solution as an initial condition. Numerical solution of (|57| ) and 
dH) at t = 5: ; WKB initial condition: — . 

of the wave function from one region to the other. The above result can also be understood 
by noting that at x = 2M the equations lose the highest order derivative responsible for 
propagation and thus the group velocity goes to zero at a; = 2M and the probability current 
also goes to zero (similar to critical layer absorption at the barrier). 

These numerical results confirm our analysis and show that when the radius of the initial 
shell is larger than its "Schwarzschild radius", the wave function is concentrated outside 
X = 2M. While when the radius of the initial shell is inside the "Schwarzschild radius" , the 
wave function leaks out very slowly. Note that the region a; ~ is always a region of small 
probability. Finally, note that these solutions show the evolution of the probability density 
of the radius of the shell measured as an intrinsic variable. As they stand, they provide no 
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FIG. 6: Numerical solution of ( |57|) and ( |5^ ) for the initial condition (^) with A = 1.0 and xq = 15.0 
for M = 5. Initial condition is outside event horizon at x = 10. Initial condition: — ; solution at 
t = 5: . Solution in < a; < 200. 
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5" 10' 15 20 ' '25'"'^"' 30 35 40 45' 50 
X 

FIG. 7: Numerical solution of ( |57|) and ( |5^ ) for the initial condition (^) with A = 1.0 and xq = 15.0 
for M = 5. Initial condition is outside event horizon at x = 10. Initial condition: — ; solution at 
t = 5: . Solution in < x < 50 showing decay behind event horizon at x = 10. 

information about the space-time geometry generated by the collapsing shell. This geometry 
must be calculated and then quantized making use of the operators generated by the motion 
of the shell. This is the subject of ongoing work. 

VI. CONCLUSIONS 

We have shown how the quantum evolution of the relativistic collapse of a shell of matter 
can be described by making use of the dimensionally reduced true relativistic Hamiltonian 
describing this process. In our analysis we have also given a procedure for quantizing a 
multi-valued Hamiltonian, and we have chosen an appropriate lapse function to produce a 
tractable Schrodinger equation. This equation has several novel features which were analyzed 
in the process of constructing the WKB solutions. These solutions have an infinite set of 
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FIG. 8: Numerical solution of ( [57| ) and (^) for the initial condition (^) with A = 1.0 and xq = 4.0 
for M = 10. Initial condition is inside event horizon at x = 20. Initial condition: — ; solution at 
t = 10: . 



bound states as well as a set of continuum eigenstates. The discrete eigenstates were shown 
to concentrate away from R = 0. Thus, in analogy to |jl5|, this suggests that the shell 
does not enter on average the singularity. The following possibility is also suggested by the 
present analysis: an incoming wave packet which classically collapses into the singularity 
R = now, due to the quantum fluctuations, is bounced back but into a new expanding 
space-time because of the small quantum fluctuations at i? ~ 0. This is a consequence of 
the Uncertainty Principle which in keeps open the throat of the wormhole and in this 
work bounces the wave function into a new expanding space-time. 

The numerical solutions show that the quantum dynamics of an initial condition, repre- 
senting a Gaussian distribution of shell radii, evolves into a sum over the discrete eigenstates. 
The critical layer nature of the process is contrasted with the behavior of other ad hoc models 
appearing in the literature which turns out to be closer to that of the hydrogen atom. 

The quantum mechanics developed in this work does not directly answer questions about 
the quantum geometry of the space. It does, however, set the framework for analyzing 
the local geometry of the space-time around the shell in terms of the operators defined in 
this work. In addition, the viewpoint we have presented provides a completely consistent 
quantization from first principles and a choice of time, related to the proper time on the 
shell, which explains in WKB terms the behavior of the dynamics observed in previous 
ad hoc models. It must be stressed that the combination of the WKB approximation and 
numerics that we used provides a complete description of a consistent quantum dynamics 
for the radius of collapsing shells of matter in a wide parameter regime. 

There are several possible directions for future work. An important study would be to 
investigate the external geometry of the system as mentioned above and the embedding of 



the shell in this spacetime. This is not an easy job. In [0, which discusses the collapse of 
null shells, it was necessary, in order to obtain a tractable problem, to define new canonical 
variables, and the results for the exterior geometry and embedding were not those that would 
seem obvious at first glance. Finding the exterior geometry would allow us to understand 
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whether the quantum mechanical expansion of the shells we found corresponds to a quantum 
mechanical tunneling through a horizon or expands into a different spacetime. Another 
direction would be to study superpositions of states with different values of M. This would 
lead to a "fuzzy" horizon. Note that in our approximation (0), these superpositions will 
have a very different evolution from the M-eigenstate we have used. 

Finally, it is important to remark that the results here obtained are within the mini- 
superspace formalism, where one first reduces the classical system and then quantizes. Re- 
cent results in isotropic (loop) quantum cosmology indicate that the inverse strategy, namely 
to first quantize end then reduce, might yield qualitative different physical predictions at 
the Planck scale [jl9|. 



APPENDIX A 



We here consider the transition layer behavior at a; = R*{E). To do this, expand Eq. 
I) about R*{E) to obtain 



or 
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l + f\R*,E)[x-R*{E)] = l + —. 



(Al) 



(A2) 



Now 9''^ is replaced by the second derivative of the boundary layer function, w, and we 
obtain 

w" + f[R*{E),E)[x- R*{E)]w. (A3) 
Since /' < 0, we have the Airy equation for w [|ll|]. It has solutions of the form 

Ai(/3i/3[x-i?*(i?)]), 



w 



/3 



- (M2 + ^2) R* + \ME'^ 



/?*2 - 2MR* 



(A4) 
(A5) 



The Airy function matches the oscillatory solution as x — R*{E) —oo, and decays ex- 
ponentially as X — R*{E) — >■ oo. In the matching of the phases of the Airy function and 
siinp2{x), we obtain the WKB energy quantization rule in the usual way . 

To consider the matching at the point R{E), we follow the same procedure. Expanding 
around R{E), we obtain 



f'[RiE)][x - RiE)]u 



1 

2n J -oo 



sinh k\e u dk 



Changing variables to x' = x — R{E) , we have 



_ 1 noo 

f'[R{E)]x'u = — I sinh fc|e*^^ u dk. 

2tT J-oo 



This equation can be solved by means of Fourier transforms to give 



. d ^ 
'dk"" 



sinh k\u. 



(A6) 



(AT) 



(A8) 
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which yields 

ii = ^e-^/o'l^i°'^«l'^, (A9) 

or 

^- |^g-.(coshfc+i)^ A;<0. ^ ^ 

Note that this function is not analytic at A; = 0. Thus, the only path of integration for 
the inversion is the real axis. We obtain 

A A 
u(x')^— e'^'' e-'^''°^'^^-^^dk + — e'^"" e'^''°^'^^~^^dk. (All) 

271 Jo 2tI J-oo 

For X > the points of stationary phase satisfy 

|sinhA;|=x', (A12) 

as they should to match the behavior for x > 0. On the other hand, for x < there are no 
points of stationary phase. The integral, after an integration by parts, shows that u decays 
as 

. . constant , /...^x 

u[x) ~ as x — > — oo. (A13) 

x' 

This is due to the nonanaliticity of | sinh k \ . 

This behavior shows us the reason for the sharp decay of amplitude observed in the 
numerics. Note that although this analysis is only local around x — R{E), it is sufficient to 
complete the WKB solution and to explain the numerics. 
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